o 


UNCLASSIFIED 


f^epAoduced 
luf.  ike 


ARMED  SERVICES  TECHNICAL  INPORMATON  AGENCY 
ARLINGTON  HALL  STATION 
AEINCTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incvirs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


ffiMORANDUM 

RM-3410-PR 


EMBER  1962 


o 


'■yi-'.».v*s!«0 

y.’r^KAirS 


CO 

<t: 


A  THIN  AIRFOIL  WITH 
A  MINIMUM  AVERAGE  COEFFICIENT  OF 
HEAT  TRANSFER  FOR  A  GIVEN  LIFT 


Translated  from  the  Russian  by  Joy  B.  Gazley 


PREPARED  FOR: 

UNITED  STATES  AIR  FORCE  PROJECT  RAND 


A  S  T  I  A 


mm 


SANTA  MONICA  •  CALIFORNIA 


MEMORANDUM 

RM-3410-PR 

NOVEMBER  1962 


A  THIN  AIRFOIL  WITH 
A  MINIMUM  AVERAGE  COEFFICIENT  OF 
HEAT  TRANSFER  FOR  A  GIVEN  LIFT 

Translated  from  the  Russian  by  Joy  B.  Gazley 


iii 


HiANSUTGR'S  PREFACE 


IhlB  Menoraadim  is  a  translation  of  tbe  Russian  paper,  "Tonkli 
Prof  11'  Krylas  Minimal 'nym  Srednlm  Koeffltslentom  Teplootdachl  Prl 
Zadannoi  Pod'*emnoi  Slle,"  (A  Thin  Airfoil  with  a  Minimum  Average  Co¬ 
efficient  of  Heat  Transfer  for  a  Given  Lift),  by  E.  V.  Bulygina,  M.  B. 
Poliakov,  and  la.  S.  Shcherbak,  published  in  Izv^stlia,  Vysshlkh 
Uchebnykh  Zavedenii,  Avlatslonnala  Tekhnika,  2,  I96I,  pp.  17-25. 

Ihls  is  one  of  a  series  of  translations  on  the  general  subject 
of  heat  transfer. 


A  THIN  AIRFOIL  WITH  A  MEMIMUM  AVERAGE  COEFFICIENT  OF  HEAT 
TRANSFER  FOR  A  GIVEN  LIFT 

E.  V.  Bulygina,  M.  B,  Poliakov,  la.  S.  Shcherfcak 

The  transition  to  high  supersonic  velocities  of  flight  leads  to  signi¬ 
ficant  kinetic  heating  of  the  wing  of  an  aircraft.  The  value  of  the  local 
heat  flow  from  the  boundary  .layer  into  the  wing  is  determined  by  the 
product  of  the  local  heat-transfer  coefficient  h^  and  the  difference  be¬ 
tween  the  local  boundary-layer  ten^jeratures  T^^  and  the  wall  tenqperatures, 

%  =  “l  M 

It  is  clear  that  a  decrease  in  the  heat-transfer  coefficient  on  the  wing 
leads  to  a  decrease  of  thermal  flows  into  the  wing.  The  coefficient  of 
heat  transfer  depends,  basically,  on  the  local  flow  parameters  at  the  outer 
edge  of  the  boundary  layer.  The  angle  of  attack  and  the  profile  shape, 
which  determine  in  addition  the  lift  coefficient  Cy,  affect  the  local  flow 
parameters. 

The  problem  set  is  to  determine  some  airfoil  shape,  which,  at  a  given 
Cy,  provides  for  a  minimum  average  heat-transfer  coefficient  at  the  profile. 
The  problem  raised  is  solved  with  the  following  sing)lifying  assumptions, 

1.  A  thin  profile  is  assumed  (the  angles  of  Inclination  of  its  upper 
and  the  lower  surface  are  the  same), 

2,  The  lower  and  the  upper  surfaces  of  the  profile  are  heat  insulated. 

5.  The  bo\mdary  layer  is  either  cort5)letely  laminar  or  con?)letely 


turbulent. 
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Let  the  equation  of  the  mean  line  of  the  profile  he: 

y  =  y(x)  (2) 

We  vrill  designate  the  angle  of  inclination  of  the  surface  of  an  infinitely 
thin  airfoil: 

e  =  y'(x) 

oince  the  profile  is  thin,  we  determine^  in  the  first  approximation,  the 
pressure  coefficient,  p,  hy  the  linear  theory: 


P  = 


'I't  -  2 


(3) 


where  is  the  Mach  nuniber  of  the  external  flow. 

The  lift  coefficient  of  the  profile  is  determined  hy  the  integral  of 
the  local  angle  of  the  mean  line  of  the  profile; 

h 

Py  =  —  r  s  (x)dx,  (4) 

where  h  is  the  chord  of  the  profile  (meters). 

The  mean  heat-transfer  coefficient  is  obtained  hy  integrating  the 
local  heat-transfer  coefficients  for  the  upper  and  lower  surfaces  of  the 
profile; 


h 


0 


where  h(x,e)  is  the  local  heat-transfer  coefficient  (k  cal/ra^  hr  °C)  and 
X  is  the  distance  along  the  chord  from  the  leading  edge  of  the  profile 
(meters).  Here,  and  below,  the  superscript  "+"  indicates  a  heat-transfer 
coefficient  for  values  of  x  and  0  on  the  lower  surface  of  the  profile  and 
the  superscript  for  values  on  the  upper  surface  of  the  profile. 

There  la  a  series  of  equations  for  the  determination  of  local  heat- 
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transfer  co-efficients  on  a  plate  tinder  a  zero  angle  of  attack. 

For  the  calculation  of  local  heat-transfer  coefficients  on  a  thin  profile 
under  a  small  angle  of  attack  it  is  possible,  in  the  first  approximation, 
to  use  the  dependencies  obtained  for  plates  under  a  zero  angle  of  attack  if 
one  substitutes  in  these  dependencies  the  local- flow  parameters  at  the 
outer  edge  of  the  profile  boundary  layer.  As  is  shown  in  (4),  this  method 
of  calculating  heat-transfer  coefficients  for  not  very  thick  profiles  with 
a  pressure  gradient  for  a  laminar  boundary  layer  gives  results  which  agree 
satisfactorily  with  calculations  by  more  accurate .methods. 

Reference  {$)  points  out  the  possibility  of  using  this  method  for  the 
turbulent  boundary  layer.  In  the  present  article  the  heat-transfer  co¬ 
efficient  is  determined  according  to  the  eq.uations  of  Ref.  (l),  lAiich, 
using  the  local  values  of  the  flow  parameters  at  the  thin  profile  in  flow, 
may  be  transformed  to: 

h  =  Afjj(H)  f^  (x)  f-  (M^,  M^),  (6) 


where 


1-n 

%  -  =PH  -h  4'“' 

(x)  = 


’h'  Pg  =  the  coefficients  of  viscosity,  temperature,  and  density  of 
air  at  an  altitude  H  at  standard  atmospheric  conditions, 

,  specific  heat  of  air  at  an  altitude  H(k  cal/kg°C)  (in 
the  calculations,  ■  0.2k), 

K  «  adiabatic  index  for  air  (  k  *  l-k), 


k 


teniperature  of  the  surface  (°K), 

\}  ■  the  local  teniperature  and  Mach  number  of  flow  at  the  outer 

edge  of  the  profile  boundary  layer. 

For  the  turbulent  boundary  layer 

A  =  14550,  n  =  0.2,  3  =  2.56,  a  =  0.44^ 
for  the  laminar  boundary  layer 

A  =  64460,  n  =  0.5,  3  =  1.9,  a  =  0.1. 

The  ten5)erature  of  the  wall  surface,  T^,  may  change  from  the  tempera¬ 
ture  of  the  external  flow,  T^,  to  the  temperature  of  the  heat-insulated 
wall; 

^a„l  •  (T) 

In  this  equation  r  is  the  tenq)erature  recovery  factor.  For  the  turbulent 
boundary  layer,  r  =  0.9,  for  the  laminar,  r  7  O.85. 

It  is  possible  to  isolate  two  extreme  cases  of  temperature  change  of 
the  wall: 


^wi  “  \wi  ^°^®spond5  to  a  continuous  steady  flight  without  taking 
into  account  heat-insulation, 


corresponds  to  instantaneous  climb  to  altitude  H  and  a 
Mach  number,  M^,  of  flight  ("instantaneous  start"). 

For  ease  of  calculation,  it  is  possible  to  combine  these  cases  if  the 
tengierature  recovery  factor  on  the  surface  is  assumed  approximately  equal 
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For  the  turbulent  boundary  layer 


7=5  for  (a  =  O.4I1), 

7  =  2.56  for  (a  =  0). 

For  the  laminar  boundary  layer 
,  .  2  for  ^ 

7  »  1.9  for  =  Tj  (a  =  0), 

The  local'-pressure  coefficient  and  local  Mach  number,  M,  with  their 
small  changes  are  combined  in  the  equation: 

r~  I  1^ 

(9) 


Mi  V(1 (l  +  |^^-i)^-l 


Thus  the  problem  of  finding  the  configuration  of  a  thin  profile  having, 
at  a  given  flight  regime,  a  minimum  meairheat- transfer  coefficient  for 
given  values  of  the  lift  coefficient  leads  to  the  finding  of  some  function 
y(x)  which  gives  the  extreme  functional  equation  (5)  for  the  given  values 
of  the  functional  equation  (4). 

For  the  solution  of  this  single  boundary-variation  problem  we  take 
advantage  of  the  LaGrange  method  of  multiples  according  to  which  the 
function  sought  must  satisfy  the  Euler  equation: 


dF  d 


dF 

\  / 


=  0 


(10) 


where 


F  =  h^(x,0)  +  h  (x,0)  +  •9(x), 


Xj^  being  the  constant  of  the  LaGrange  multiple.  Since  F  apparently  does 
not  depend  on  y,  Euler's  equation  is  reduced  and  then  its  integral  is 
found; 


■ 

df-  dM,  d- 

+ 

1  df-  dM,  d-, ' 

" 

P  .  i  Pi 

+ 

p  .  i  .  Pi 

5m:  5^  35" 

1  Pi 

i  i  Pi 

=  C 


(11) 
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The  values  ^  ,  which  are  calculated  on  the  basis  of  Eq,  (8),  enter 
Eq.  (ll): 


On  the  basis  of  Eq.  (9)  we  obtsdn 


The  derivative  ^  for  the  upper  surface  will  have  the  value; 


for  the  lower  surface: 


The  relations 


(I4a) 


(14b) 


for  the  upper  and  lower  surfaces  have  one  and  the  seune  form  (12),  (13) 
they  differ  only  in  the  fact  that  in  the  first  case  the  pressvire  coefficient 
is  calculated  as  positive,  in  the  second  as  negative. 

The  equation  which  combines  the  pressure  p  and  the  coordinate  of 


point  X  takes  the  form: 
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-a 


l-2x 


(1-n)  +  (x-1)  -  7) 


i»x 


(1  +  ^^^)  (1  +  |f^  i)  ^  -  1  + 


7+1 


l-2x 

(1  +  I  5)  ^  (1-n)  +  (x-1)  (^  -  7) 


(1  +  ^Rg)  (1  +  5)  ""  -  1 


(1  +  ^M?) 


7+1 


(15) 


Here  the  values  n,  7  are  the  same  as  in  Eq.  (8), 


=1  = 


v  €  •  ^ 


V2(x-1)  .  Afj  (1  t 


@+l 


This  equation  defines  p  as  an  implicit  function  of  x: 


0(p)  =  0(' 


26  \  .  n.n 


V'<  -  1 


where  K  =  C^n  is  the  new  derivative  constant. 


(16) 


It  isn't  difficult  to  see  from  Eq.  (I5),  for  p  =  0  (or  6  =  O),  that 
the  function  0  (p)  reduces  to  zero  and,  consequently,  for  any  value  of  the 
derivative  constant  C^,  the  angle  of  inclination,  equal  to  zero,  corresponds 
to  the  origin  of  the  coordinates.  Eq.  (I5)  permits  the  easy  determination 
of  X  as  a  function  of  0. 
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In  the  nature  of  an  example ,  In  Figs.  1  and  2  axe  Introduced  graphs 
of  this  dependency  calculated  for  conditions  of  "Instantaneous  start"  and 
subsequent  flight  at  an  altitude  H  =  20  km  with  a  Mach  number  =  3  for 
the  laminar  and  the  turbulent  boundary  layers. 

Hie  lift  coefficient  Cy  is  obtained  by  numerous  integrations  of  the 
curve  p  (x)  (or  0(x)): 


The  coordinate  of  the  profile  y  is  determined  by  the  integration  of 

the  angle  of  inclination  of  the  mean  line  of  the  profile: 
xk 

^y  =  J  0(xX)  d(x\)  . 

The  lift  coefficient  Cy  and  the  final  ordinate  of  the  profile  y^^^  are 
combined  by  the  simple  dependency: 


k 


A  specific  value  of  the  chord  b\  corresponds  to  each  value  of  Cy. 

The  dependence  of  the  rate  b\  on  Cy  for  =  3  and  H  =  20  km  is 
represented  for  the  turbulent  and  the  laminar  boundary  layers  in  Figs.  3 
and  h. 

The  curves  of  the  corresponding  profiles  at  a  given  value  of  Cy  =  O.178 
are  shown  in  Fig.  3> 

As  is  evident  from  this  figure,  the  profiles  essentially  differ  from 
the  plate  to  the  adjacent  leading  surface  where  the  profile  obtained  has  a 


9 


zero  angle  of  inclination. 

It  niust  be  noted  that  the  value  of  the  heat-transfer  coefficient,  as 
in  the  case  of  the  turbulent,  so  also  in  the  case  of  the  laminar  boundary 
layer,  depends,  essentially,  on  the  temperature  factor. 

But  since  the  form  of  the  profile  for  the  case  of  tvo  edges  follows 
from  the  calculations  and  from  Fig.  5,  and  are  obtained 

as  the  same  tiling;  i.e.,  the  form  of  the  optimal  profile  depends  slightly 
on  the  teii5>erature  factor.  Thus,  the  profile  found  will  be  optimal  also 
at  all  intermediate  wall  ten^jeratures. 


Pig.  2.  Die  Lftcilnar  Boundary  Layer 
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